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The nuclear mass dependence of the number of short-range correlated (SRC) proton-proton (pp) 
and proton-neutron (pn) pairs in nuclei is a sensitive probe of the dynamics of short-range pairs in 
the ground state of atomic nuclei. This work presents an analysis of electroinduced single-proton 
and two-proton knockout measurements off ^^Al, ®®Fe, and ^°®Pb in kinematics dominated 
by scattering off SRC pairs. The nuclear mass dependence of the observed A{e,e'pp)/^^C{e,e'pp) 
cross-section ratios and the extracted number of pp- and pn-SRC pairs are much softer than the 
mass dependence of the total number of possible pairs. This is in agreement with a physical picture 
of SRC affecting predominantly nucleon-nucleon pairs in a nodeless relative-S state of the mean-field 
basis. 

PACS numbers: 25.30.Rw, 25.30.Fj, 24/10.-i 


Introduction: The nuclear momentum distribution 
(NMD) is often quoted as being composed of two 
separate parts m- Below the Fermi momentum 
{kp ^ 250 MeV/c) single nucleons move as independent 
particles in a mean field created by their mutual inter¬ 
actions. Above the Fermi momentum (fc > kp) nucleons 
predominantly belong to short-range correlated (SRC) 
pairs with high relative and low center-of-mass (c.m.) 
momenta, where high and low are relative to the Fermi 
momentum mi]. In addition to its intrinsic interest, 
the NMD and its division into mean-field and correlated 
parts is relevant to two-component Fermi systems [S], 
neutrino physics [B E], and the symmetry energy of 
nuclear matter [12]. 

The mean-held and long-range aspects of nuclear dy¬ 
namics have been studied extensively since the dawn of 
nuclear physics. The effect of long-range correlations on 
the NMDs is limited to momenta which do not extend 
far beyond kp [13]. Study of the short-range aspects of 
nuclear dynamics has blossomed with the growing avail¬ 
ability of high-energy high-intensity electron and pro¬ 
ton accelerators. Recent experiments conhrm the pre¬ 
dictions that SRC pairs dominate the high-momentum 
tails (fc > kp) of the NMDs [IHZ], accounting for 20-25 
% of the NMD probability density [MUTT] . These high- 
momentum tails have approximately the same shape for 
all nuclei [21 [31 [HI ICTITB] , differing only by scale factors 
which can be interpreted as a measure of the relative 
number of SRC pairs in the different nuclei. In this work, 
we aim at understanding the underlying dynamics which 
give rise to this universal behavior of the high-momentum 
tail. 

An intuitive picture describing the dynamics of nuclei 
including SRCs is that of independent bound nucleons 
moving in the nucleus, occasionally getting sufficiently 
close to each other to temporarily fluctuate into SRC- 
induced nucleon-nucleon pairs. This picture can be for¬ 


mally implemented in a framework in which one shifts the 
complexity of the nuclear SRC from the wave functions to 
the operators by calculating independent-particle model 
(IPM) Slater determinant wave functions and acting on 
them with correlation operators to include the effect of 
SRCs [TBfl20j . The observed number of proton-proton 
(pp) and proton-neutron (pn) SRC pairs in various nuclei 
can then be used to constrain the amount and the quan¬ 
tum numbers of the initial-state IPM nucleon-nucleon 
(SRC-prone) pairs that can fluctuate dynamically into 
SRC pairs through the action of correlation operators. 

In this paper, we will extract the relative number of 
pp-SRC and pn-SRC pairs in different nuclei from mea¬ 
surements of electroinduced two-proton and one-proton 
knockout cross-section ratios for medium and heavy nu¬ 
clei (^^Al, ®®Fe, and ^°®Pb) relative to ^^C in kinematics 
dominated by scattering off SRC pairs mm- In these 
kinematics in the plane-wave approximation A(e, e'pp) 
cross sections are proportional to the number of pp-pairs 
in the nucleus and A(e, e'p) cross sections are propor¬ 
tional to twice the number of pp pairs plus the number of 
pn pairs (2pp-|-pn). Therefore, after correcting the mea¬ 
sured cross sections for rescattering of the outgoing nu¬ 
cleons from the residual nucleus (final state interactions 
or FSI), the relative number of pp and pn pairs will be 
extracted from measurements of A{e,e'pp)/^'^C{e,e'pp) 
and A(e, e'p)/^^C(e, e'p) cross-section ratios [5]. 

We will then compare the A{e, e'pp)/^^C(e, e'pp) cross- 
section ratios and the extracted number of pp and pn 
pairs to factorized calculations using different models of 
nucleon pairs in order to deduce the quantum numbers 
of the IPM SRC-prone pairs. We will provide strong evi¬ 
dence that the relative quantum numbers of the majority 
of the SRC-susceptible pairs are ^So(l) for pp and ^Si(O) 
for pn. Hereby, we used the notation ^"'+^Ls(T) to iden¬ 
tify the pair’s quantum state (T is the total isospin). 

Experiment: The one- and two-proton knockout mea- 
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FIG. 1: (color online). The distribution of the cosine of the 
angle between the missing momentum of the leading proton 
and the recoil proton for The insert shows the same dis¬ 
tribution for ^^Al, ®®Fe, and ^°®Pb. The dashed (red on¬ 
line) line shows the distribution of the random phase-space 
extracted from mixed-events. 


surements analyzed in this paper were described in [S] 
and its supplemental information. They were carried 
out using the CEBAF Large Acceptance Spectrometer 
(CLAS) [12], located in Hall-B of the Thomas Jefferson 
National Accelerator Facility (Jefferson Lab) in Newport 
News, Virginia. The data were collected in 2004 us¬ 
ing a 5.014 GeV electron beam incident on ^^Al, 
®®Fe and ^°®Pb targets. The scattered electron and 
knocked out proton(s) were measured with CLAS. We 
selected A(e, e'p) events in which the electron interacts 
with a single fast proton from a SRC nucleon-nucleon 
pair in the nucleus by requiring large four-momentum 
transfer {Q^ > 1.5 GeV^), Bjorken scaling parameter 

= 2 mHU} ^ missing momentum 300 < 

IPmissI < 600 MeV/c. The four-momentum transfer 
= q-q — (^) where q and w are the three-momentum 
and energy transferred to the nucleus respectively; to at is 
the nucleon mass; the missing momentum Pmiss = Tp — (f, 
and Pp is the knockout proton three-momentum. We also 
required that the knockout proton was detected within 
a cone of 25° of the momentum transfer q and that it 
carried at least 60% of its momentum (i.e. > 0.6). 

To suppress contributions from inelastic excitations of 
the struck nucleon we limited the reconstructed missing 
mass of the two-nucleon system TOmiss <1.1 GeV/c^. 

The A(e, e'pp) event sample contains all A{e, e'p) 
events in which a second, recoil, proton was detected 
with momentum greater than 350 MeV/c. Fig. shows 
the distribution of the cosine of the angle between the 
initial momentum of the knockout proton and the recoil 
proton for these events |S]. The recoil proton is emitted 
almost diametrically opposite to the missing-momentum 


direction. The observed backward-peaked angular distri¬ 
butions are very similar for all nuclei and are not due to 
acceptance effects as shown by the angular distribution 
of mixed events. These distributions are a signature of 
scattering on a nucleon in a SRG pair, indicating that 
the two emitted protons were largely back-to-back in the 
initial state, having large relative momentum and small 
c.m. momentum |6l|23|. Further evidence of scattering on 
a SRG nucleon pair is that the recoil proton was emitted 
at forward angles (i.e., angles in the range 20-60° with 
respect to q). 

The A{e,e'p)/^^C{e,e'p) and A{e,e'pp)/^'^C{e,e'pp) 
cross-section ratios are obtained from the ratio of the 
measured number of events, normalized by the incident 
integrated electron flux and the nuclear density of each 
target. During the experiment all solid targets were 
held in the same location, the detector instantaneous 
rate was kept constant, and the kinematics of the mea¬ 
sured events from all target nuclei were almost identi¬ 
cal mm- Therefore detector acceptance effects cancel 
almost entirely in the A{e,e'pp)/C{e,e'pp) cross section 
ratios. Due to the large acceptance of GLAS, radiative 
effects affect mainly the electron kinematics. These cor¬ 
rections were calculated in Ref. m for the extraction of 
the A(e, e'p)fG{e, e'p) cross section ratio. As the electron 
kinematics is the same for the A(e, e'p) and A(e, e'pp) re¬ 
actions, the same corrections are used here to extract the 
A{e, e'pp)/C{e, e'pp) cross-section ratios. See Ref. [8| for 
additional details. 

FSI model: In order to extract the underlying relative 
number of pp and pn SRG pairs in nuclei from the mea¬ 
sured cross-section ratios, we must correct the data for 
FSI effects |8]. Alternatively, in order to compare the 
measured ratios to calculations, we must correct either 
the data or the calculation for FSI effects. The two dom¬ 
inant contributions are: (1) attenuation of the outgoing 
nucleon(s) upon traversing the residual A — 1 or A — 2 
nucleus, and (2) rescattering of a neutron into a proton 
(i.e., single charge-exchange (SCX)). SGX can lead to a 
pp final state which originates from a pn pair. 

The effect of FSIs of the ejected pair with the re¬ 
maining A — 2 spectators was computed in a relativis¬ 
tic multiple-scattering Glauber approximation (RMSGA) 
pn E5] . The RMSGA is a multiple-scattering for¬ 
malism based on the eikonal approximation with spin- 
independent NN interactions. We have included both the 
elastic and the SGX rescattering of the outgoing nucleons 
with the A —2 spectators. The three parameters entering 
in the RMSGA model are taken from NN scattering data 
and yield an excellent description of the world’s A(e, e'p) 
transparency data [IS]. The RMSGA caclulations yield 
attenuation coefficients that are similar to the power-law 
results obtained in analysis of nuclear transparency mea¬ 
surements m- In this work no free parameters are tuned 
to model the FSI effects in the A(e, e'p) and A(e, e'pp) 
data under study. 
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The sex probabilities are calculated in a semi-classical 
approximation. The probability of charge-exchange re¬ 
scattering for a nucleon with initial IPM quantum num¬ 
bers a which is brought in a continuum state at the co¬ 
ordinate r is modeled by, 


pa(/3) 

^CX 


(r ) 


^- 1-00 

= 1 - exp[-crcx(s) / Az'p'^^iz')]. 

J Z 


( 1 ) 


The z-axis is chosen along the direction of propagation of 
the nucleon with initial quantum numbers a. The quan¬ 
tum numbers of the correlated partner in the SRC pair 
are denoted with /?. The is the density of the resid¬ 
ual nucleus available for SCX reactions. Obviously, for an 
ejected proton (neutron) only the neutron (proton) den¬ 
sity of the residual nucleus affects SCX reactions, crcx(s) 
in Eq. Q, with s the total c.m. energy squared of the 
two nucleons involved in the SCX |26], can be extracted 
from elastic proton-neutron scattering data m- 

Cross-section Model: As outlined in Refs. [231128], in 
the spectator approximation it is possible to factorize 
the A(e, e'pN) cross section in kinematics probing short- 
range correlated pairs as 


d®cr [A(e, e'pN)] 
Pi2(i^ki2 


KepNaepN{h2)Ff^^\Pl2), ( 2 ) 


where fie' is the solid angle of the scattered electron, 
and fci 2 and P 12 are the relative and c.m. momenta of 
the nucleon pair that absorbed the virtual-photon. The 
KepN is a kinematic factor and aepN{ki 2 ) is the cross 
section for virtual-photon absorption on a correlated pN 
pair. The F^^^\Pi 2 ) is the distorted two-body c.m. 
momentum distribution of the correlated pN pair. In the 
limit of vanishing FSIs, it is the conditional c.m. momen¬ 
tum distribution of a pN pair with relative S„=o quan¬ 
tum numbers. Distortions of F^^^\Pi 2 ) due to ESI are 
calculated in the RMSGA. The factorized cross-section 
expression of Eq. ([^ hinges on the validity of the zero- 
range approximation (ZRA), which amounts to putting 
the relative pair coordinate ri 2 to zero. The ZRA works 
as a projection operator for selecting the very short-range 
components of the IPM relative pair wave functions. 

The probability for charge-exchange reactions in pN 
knockout is calculated on an event per event basis, us¬ 
ing the SRC pair probability density F^^^\Ri 2 ) in the 
ZRA corrected for ESI. With the aid of the factorized 
cross-section expression of Eq. the phase-space inte¬ 
grated A(e, e'pN) to ^^C(e, e'pN) cross-section ratios can 
be approximately expressed as integrals over distorted 
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FIG. 2: (color online). The mass dependence of the 
A{e, e'pp)/^'^C{e, e'pp) cross-section ratios. The points show 
the measured, uncorrected, cross section ratios. The lower 
orange band and upper grey line denote ZRA reaction-model 
calculations for ^^Al, ®®Fe, and based on Eq. (j^ 

with and without FSI corrections respectively. The width of 
the ZRA-RMSGA band reflects the maximum possible effect 
of SGX. 


c.m. momentum distributions, 

cr [A(e, e'pN)] 
a [^^C(e, e'pN)] 

f d^ne>d^ki2N,p^aepN(kl2) f d^Pl2F^^^^\Pi2) 

f d^n,,d^ki2N,^Na,pM(ki2) f d^Pl2F^^^^\Pi2) 

f d^Pl2Ff(^\P,2) 
f d^Pl2F^^^^\P,2) ' 

In the absence of FSI, the integrated c.m. momentum 
distributions / d^Pi 2 F^^^^\Pi 2 ) equal the total num¬ 
ber of SRC-prone pN pairs in the nucleus A. Hence, 
the cross section ratios of Eq. (§ provide access to the 
relative number of SRC pN-pairs up to corrections stem¬ 
ming from FSI. We have evaluated the ratios of the dis¬ 
torted c.m. momentum distributions of Eq. ([^ over 
the phase space covered in the experiment. Given the 
almost 47r phase space and the high computational re¬ 
quirement of multidimensional FSI calculations, we use 
an importance-sampling approach. The major effect on 
the c.m. momentum distribution F^^^^\Pi 2 ) when in¬ 
cluding FSIs is an overall attenuation, the shape is almost 
unaffected [53]. Motivated by this, we used the c.m. mo¬ 
mentum distributions without FSI as sampling distribu¬ 
tion for the importance sampling in the FSI calculations. 
When convergence is reached, the computed impact of 
FSI is extrapolated to the whole phase space. 

Results: Figure shows the measured uncorrected 

ap^c^eei^p)] ccoss-section ratios compared with the ZRA 
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reaction-model calculation with and without RMSGA 
FSI corrections. The first striking observation is that the 
measured cross-section ratios increase very slowly with A 
(e.g., the Pb/C ratio is only 3.8±0.5). For contrast, com¬ 
binatorial scaling based on the number of pp pairs leads 
to a ratio of over 200. The ZRA-RMSGA calculations 
agree well with the measured data, yielding a Pb/G ratio 
of 4.96lg :[ 4 . The ZRA and ZRA-RMSGA calculations 
assume that only pairs with a finite probability density 
at relative coordinate zero contribute to the cross-section. 
This is consistent with assuming that only IPM pairs in 
a relative S„=o state contribute. 

Figure j^shows the number of pp- and pn-SRG pairs in 
various nuclei relative to Garbon extracted from the mea¬ 
sured A{e,e'pp)/C{e,e'pp) and A{e,e'p)/C{e,e'p) cross- 
section ratios following the method outlined in Ref. [5] 
with RMSGA corrections for FSI and SGX. The ex¬ 
tracted number of pp pairs are very sensitive to SGX. 
If the virtual photon is absorbed on a pn pair and the 
neutron subsequently undergoes a single charge exchange 
reaction with a proton, two protons will be detected in 
the final state. These events must be subtracted in order 
to extract the number of pp-SRC pairs. As the contri¬ 
bution from these pn pairs to the pp final state is com¬ 
parable to the number of initial pp pairs, this leads to 
a large uncertainty in the number of pp pairs, especially 
for heavy nuclei. 

Figure [^also shows the expected number of pp and pn 
SRG pairs relative to Garbon for different quantum num¬ 
bers of the IPM pairs that can dynamically form SRG 
pairs through the action of correlation operators. These 
include (a) all possible NN pairs (i.e. Z(Z-l)/(6 • 5) and 
ZN/(6 • 6) for pp and pn pairs respectively), (b) pairs in a 
nodeless relative S state (i.e. S„=o)j and (c) L < 1 pairs 
(i.e. both S and P state pairs). Those ”S„=o” pairs are 
characterized by the (n = 0, L = 0) quantum numbers for 
their relative orbital motion. Of all possible states for the 
pairs, the S„=o pairs have the highest probability for the 
two nucleons in the pair to approach each other closely. 
Close-proximity IPM pn pairs in a ^Si(O) state are highly 
susceptible to the tensor correlation operator that cre¬ 
ates SRC pairs in a spin-triplet state with predominantly 
deuteron-like quantum numbers (L = 0, 2; T = 0; S' = 1). 

We determine the number of pairs in each case us¬ 
ing an IPM harmonic-oscillator basis and performing a 
standard transformation to relative and center-of-mass 
coordinates as detailed in Ref. [29] . The relative number 
of pairs are displayed in Fig. [^ and listed in Table [^ As 
can be seen, both (a) the naive combinatorial assump¬ 
tion and (c) the calculations that include IPM S and P 
pairs contributions both drastically overestimate the in¬ 
crease in the number of pairs with A. The ZRA and S„=o 
pairs counting calculations are in fair agreement with the 
extracted number of pp and pn pairs. 

As both the ZRA and the S„=o pair counting project 
IPM states onto close-range pairs, we expect the two 
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TABLE I: The relative number of SRC pp and pn pairs calcu¬ 
lated using Sn=o counting and the ZRA reaction model com¬ 
pared to the extracted values from the measured A{e, e'p) and 
A(e, e'pp) ratios after correcting for FSI effects. The error in¬ 
cludes the uncertainties on the cross-section ratios and FSI 
calculations. 

methods to produce a similar mass dependence of the 
number of SRG pairs The ZRA predicts a somewhat 
softer mass dependence (oc vs 

This can be explained by the fact that the ZRA is a 
more restrictive projection on close-proximity pairs than 
the S„=o counting which accounts also for ri 2 yf 0 con¬ 
tributions. 

The observed agreement with the experimental data 
indicates that correlation operators acting on IPM S„=o 
pairs are responsible for the largest fraction of the high- 
momentum nucleons in nuclei. This gives further support 
to the assumption that the number of IPM pairs with 
quantum numbers S„=o is a good proxy for the number 
of correlated pairs in any nucleus A HHUminOI. This is 
also consistent with an analysis of the cross section of the 
ground-state to ground-state transition in high-resolution 
^®0(e, e'ppy^C measurements |3T1|32| which provided ev¬ 
idence for the ^So(l) dominance in SRG-prone pp pairs. 

Conclusions: We have extracted the relative num¬ 
ber of np and pp SRG correlated pairs in nucleus A 
relative to Carbon from previously published measured 
A(e, e'pp)fC(e, e'pp) and A(e, e'p)/C(e, e'p) cross section 
ratios corrected for final state interactions. The relative 
number of np and pp pairs increases much more slowly 
with A than expected from simple combinatorics. 

We calculated the cross section in a framework which 
shifts the complexity of the nuclear SRC from the wave 
functions to the operators by calculating independent- 
particle model (IPM) Slater determinant wave functions 
and acting on them with correlation operators to include 
the effect of SRCs [TStiSU] . This framework is also well 
suited for explaining the recently observed similarities 
between the contact term measured in ultra-cold atomic 
systems and that extracted for atomic nuclei [9|. The 
uncorrected A{e,e'pp)/C{e,e'pp) cross section ratios are 
consistent with a zero range approximation (ZRA) cal¬ 
culation including the effects of FSI. 

Due to factorization, the ratio of calculated cross sec¬ 
tions is approximately equal to the ratio of the distorted 
c.m. momentum distributions. In the absence of FSI, the 
integrated c.m. momentum distribution equals the total 
number of SRC-prone pairs in that nucleus. We com- 
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FIG. 3: (color online). The mass dependence of the num¬ 
ber of pp (top panel) and pn (bottom panel) SRC pairs of 
nucleus A relative to Data are extracted from the mea¬ 
sured CLAS A{e, e'p) and A{e, e'pp) cross section ratios [51l21| 
after correcting for FSI. Error bars include the estimated un¬ 
certainty on the cross-section ratios and the FSI corrections. 
The green sqnares correspond with unconditional counting of 
the pp pairs i.e. (Z(Z-l)/30 in the upper panel) and pn pairs 
(ZN/36 in the bottom panel) for the nuclei ^^Al, 

^°Ca, ^®Ca, ®®Fe, ®®Cu, i°®Ag and ^“^Pb. The yellow dia¬ 
monds are the ratios obtained by counting IPM pairs in a rel¬ 
ative S and P state. The blue triangles count IPM Sn=o pairs. 
The solid line denotes the result of a reaction-model calcula¬ 
tion for scattering from close-proximity pairs (Eq. <§) which 
takes full account of the experimental phase space. This cal¬ 
culation does not include FSI corrections as these are applied 
to the data, see text for details. 


pared three choices of SRC-prone pairs to the data: (a) 
all pairs, (b) pairs in a nodeless relative S state (S„=o)j 
and (c) L < 1 pairs (i.e., both S and P). 

We found that the soft mass dependence of the mea¬ 
sured A{e, e'pp) cross-section ratios agrees with scatter¬ 
ing from highly selective close-proximity pairs (i.e., only 
IPM relative S„=o pairs). The mass dependence of the 
extracted ratios of the number of short-range correlated 
pp and pn pairs provides additional support for this con¬ 
clusion. All these results consistently hint at a physical 
picture whereby the aggregated effect of SRC in the nu¬ 
clear wave function is determined to a large extent by 
mass-independent correlation operators on S„=o pairs. 
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